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Abstract 


The four-lepton { = e,/j.) production cross section is measured in the mass range from 
80 to 1 000 GeV using 20.3 fb“* of data in pp collisions at = 8 TeV collected with the 
ATLAS detector at the LHC. The events are produced in the decays of resonant Z and 
Higgs bosons and the non-resonant ZZ continuum originating from qq, gg, and qg initial 
states. A total of 476 signal candidate events are observed with a background expectation of 
26.2+3.6 events, enabling the measurement of the integrated cross section and the differential 
cross section as a function of the invariant mass and transverse momentum of the four-lepton 
system. 

In the mass range above 180 GeV, assuming the theoretical constraint on the qq production 
cross section calculated with perturbative NNLO QCD and NLO electroweak corrections, 
tbe signal strength of the gluon-fusion component relative to its leading-order prediction is 
determined to be pgg = 2.4 +1.0 (stat.) + 0.5 (syst.) + 0.8 (theory). 
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Reproduction of this article or parts of it is allowed as specihed in the CC-BY-3.0 license. 




1 Introduction 


This paper presents measurements of the production of four isolated charged-leptons in proton-proton 
collisions at a centre-of-mass energy of -^ = 8 TeV using 20.3 fb“^ of data collected with the ATLAS 
detector at the LHC. For the four-lepton {A£, i = e,fi) production, both the integrated cross section and 
the differential cross sections as functions of invariant mass {m 4 () and transverse momentum {pj) of 
the system, are measured in a mass range 80 < < 1 000 GeV. In addition, the signal strength 

of gluon fusion {ggF) production relative to its leading-order (LO) QCD estimate is measured. These 
measurements test the validity of the Standard Model (SM) through the interplay of QCD and electroweak 
effects for different At production mechanisms as described by the LO Feynman diagrams shown in 
Figure 1. 



Figure 1; The LO Feynman diagrams for the qq- and grg-initiated production of A€\ (a) i-channef production of 
qq FF with associated radiative decays to an additionaf fepton pair; (b) f-channei production of 

qq ^ A(\ (c) Higgs-boson production through gfuon fusion gg —> H^*'> —> ZZ*** —> A{\ (d) non-resonant 

At production through the quark-box diagram gg z**^Z*** ^ At. The notation stands for production of on- 
and off-sheff Z bosons (Z and Z*) and production of off-sheff photons (y*). 

The At signal events come from the decays of resonant Z and Higgs bosons and the non-resonant ZZ 
continuum produced from qq, gg, and qg initial states, which are briefly discussed below. 

• ^^-initiated At production 

The tree-level diagrams for qq At production are shown in Figure 1(a) and Figure 1(b). The cross 
section as a function of is shown in Figure 2 (the dashed black histogram). The At event production 
at the Z resonance occurs predominantly via the ^-channel diagram as shown in Figure 1(a), and was 
measured previously by the ATLAS and CMS collaborations [1, 2]. In the At invariant mass region above 
the Z resonance the At event production mainly proceeds through the t-channel process as shown in 
Figure 1(b). The cross section significantly increases when both Z bosons are produced on-shell, resulting 
in a rise in the mm spectrum around 180 GeV. In addition, a small portion of the At events with the qq 
initial state can be produced from the vector-boson scattering {VBS) process. 

• ^gr-initiated At production 

The LO diagrams of the Higgs-boson production and non-resonant At production via ggF are shown in 
Figure 1(c) and Figure 1(d), respectively. The cross sections as a function of mm are shown in Figure 2 
(the colored histograms). The features of the At events from the decays of Higgs-boson and continuum 
ZZ production via ggF are described below. 

(1) The dominant Higgs-boson production mechanism is ggF. Other Higgs-boson production mechan¬ 
isms, vector-boson fusion (VBF), vector-boson associated production (VH), and top-pair associated pro- 
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Figure 2: The differential cross sections, dcrldnim versus the invariant mass of the four leptons m^(, calculated by 
MCFM from the qq and gg initial states at yfs-% TeV for the lelji hnal state in the experimental fiducial phase 
space (see Table 2 for definition). The inclusive gg —> Al distribution is the sum of the gg ^ H and the 

gg —» ZZ ^ 4{, and interference terms. The calculation of the qq —> differential production cross section 
includes perturbative QCD corrections at NLO, while the distributions from the gg initial state are calculated at LO. 
The NNLO ^T-factors are applied to on-shell Higgs-boson production. 


duction {ttH), contribute less than 15% to the on-shell Higgs-boson decay to ZZ* event rate. The on-shell 
Higgs-boson production and decay leads to a narrow resonance around 125 GeV, which has been a key 
signature in the Higgs-boson discovery by the ATLAS [3] and CMS [4] collaborations. The off-shell 
Higgs-boson production has a large destructive interference with continuum ZZ production from the ggF 
processes [5-7]. This effect can be observed in the high-mass tail of the distributions shown in Figure 2, 
and has been used as a tool to constrain the total Higgs-boson width by the ATLAS and CMS collabora¬ 
tions [8, 9]. 

(2) The non-resonant ZZ —> A£ production via ggF, includes the production of off-shell Higgs bosons and 
continuum ZZ production as well as their interference. This process produces a sizeable number of Ai 
events in the > 2xmz mass region and dominates the total ^g-initiated 4£ production. 

Contributions from different processes have different shengths as a function of (Figure 2) and 
Therefore, differential A{ production cross sections are measured separately as a function of and 
p^. The measurement of the integrated cross section is first performed in the experimental fiducial phase 
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space, and then extended to a common phase space for three channels: Ae, 4/r, and 2e2/r. This common 
phase space is defined by 80 < < 1 000 GeV, me+i- > 4 GeV, pj.' >2 GeV, and the presence of four 

leptons each with pr > 5 GeV and \r]\ < 2.8. 

Currently, the gluon-fusion production is estimated theoretically with only a LO QCD approximation for 
the gg continuum production [ 6 , 10]. In this analysis the mass range above 180 GeV is used to determine 
the signal strength of the gluon-fusion component with respect to its LO prediction. This is done by fitting 
the observed m 4 ( spectrum using the next-to-next-to-leading-order (NNLO) QCD theoretical prediction, 
corrected for next-to-leading-order (NLO) electroweak effects, for the production originating from the qq 
initial state. 


2 The ATLAS detector 

The ATLAS detector [11] has a cylindrical geometry^ and consists of an inner tracking detector (ID) 
surrounded by a 2 T superconducting solenoid, electromagnetic and hadronic calorimeters, and a muon 
spectrometer (MS) with a toroidal magnetic field. The ID provides fracking for charged particles for 
\q\ < 2.5. If consisfs of silicon pixel and sfrip defectors surrounded by a sfraw fube fracker fhaf also 
provides fransifion radiafion measuremenfs for election idenfificafion. The electiomagnefic and hadronic 
calorimeter system covers fhe pseudorapidify range \q\ < 4.9. For \q\ < 2.5, fhe liquid-argon elecfromag- 
nefic calorimefer is finely segmented and plays an imporfanf role in fhe election idenfificafion. The MS 
includes fasl tiigger chambers {\q\ < 2.4) and high-precision fracking chambers covering \q\ < 2.1. A 
fhree-level tiigger system selecfs evenfs fo be recorded for offline physics analysis. 


3 Signal and background simulation 

The signal modelling for qq AC producfion uses fhe POWHEG-BOX Monfe Carlo (MC) program [12- 
14], which includes perfurbafive QCD corrections af NLO. The producfion fhrough fhe qg inifial sfafe is an 
NLO confribufion fo fhe qq process. The CTIONLO [15] sef of parfon disfribufion functions (PDFs), wifh 
QCD renormalisafion and faclorisafion scales (/ir, /ip) sef to are used to calculate fhe cross section 
and generafe fhe kinematic disfribufions. The NNLO QCD [16] and fhe NLO electioweak (EW) [17] 
corrections are applied to fhe NEO cross section calculated by POWHEG-BOX as a function of fhe AC 
mass for fhe kinematic region where bofh Z bosons are produced on-shell. Following fhe same approach 
as described in Ref. [ 8 ], fhe AC evenf disfribufions are re-weighfed fo mafch fhose expected when using 
QCD scales of m 4 ^/ 2 . This is done to unify fhe QCD scales used in simulafion of fhe qq and fhe gg 
processes. 

The signal modelling of fhe on-shell Higgs-boson producfion via fhe ggF and VBF mechanisms uses 
POWHEG-BOX which provides calculations af NEO QCD, wifh fhe CTIONEO PDFs and pK,pv = " 14 ^. 
The Higgs-boson producfion via fhe VFI and tlH mechanisms is simulated wifh PYTHIA 8 [18]. The 
NNEO QCD and NEO EW effecls on fhe cross-secfion calculations for on-shell Higgs-boson production 


* ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector 
and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the y-axis points 
upward. Cylindrical coordinates (r, (p) are used in the transverse plane, ^ being the azimuthal angle around the beam pipe. 
The pseudorapidity is defined in terms of the polar angle 6 as - In tan(6/2). 
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are summarised in Ref. [19]. The expected event yields of on-shell Higgs boson are normalised to the 
higher-order corrected cross sections. 

The non-resonant signal production includes off-shell Higgs-boson production, continuum ZZ pro¬ 
duction, and their interference. The LO MCFM generator [20] is used to simulate the non-resonant ggF 
production, with the CTIONNLO [21] set of PDFs with QCD scales of /rR,/rF set to while the 

LO MADGRAPH generator [22] is used to simulate non-resonant VBF and VBS production and their 
interference. The NNLO QCD corrections are available for off-shell Higgs-boson production [23] and 
for the interference between off-shell Higgs bosons and ZZ pairs from the gg initial state [24]. However, 
no higher-order corrections are available for the continuum gg ZZ process, which dominates the 
events from the gg initial state in the region outside the Higgs-boson resonance. Therefore, the LO cross 
section is used for the normalisation of the events produced in gluon-fusion processes. 

All the signal MC generators are interfaced to PYTHIAS for parton shower simulation, except MADGRAPH, 
which is interfaced to PYTHIA6 [25]. 

Backgrounds in this analysis include reconstructed A€ events from Z-rjets, tt, diboson {ZW, Zy and double 
Drell-Yan), triboson VW {V=Z, W), and VFI (H W), and Z-i-top (tt and t) processes, which are also 
simulated. 

The reducible background fromZ-i-jets production, which includes light- and heavy-flavour contributions, 
is modelled using both SHERPA [26] and ALPGEN [27]. The Zy process is simulated with SHERPA. 
The ft background is modelled using POWHEG-BOX. 

Background events from ZH production, where Z ^ and FI —> VV (W = WW or ZZ with two leptons 
and two neutrinos or two leptons and two jets in the final state), are simulated with PYTHIAS. The ZW 
and the tZ processes are simulated with SHERPA and MADGRAPH, respectively. The irreducible back¬ 
ground from WV and ttZ is modelled with MADGRAPH. Einally, the double-Drell-Yan ZZ production 
is modelled with PYTHIAS. 

Eor background modelling the POWHEG-BOX and MADGRAPH generators are interfaced to PYTHIA6 
for the parton shower, hadronisation and underlying-event simulation. The AEPGEN generator is inter¬ 
faced to HERWIG [2S] for the parton shower and to JIMMY [29] for the underlying event simulation. 
SHERPA uses built-in models for both the parton shower and underlying-event description. 

Both the signal and background MC events are simulated using the ATEAS detector simulation [30] based 
on the GEANT4 [31] framework. Additional pp interactions in the same and nearby bunch crossings 
(pile-up) are included in the simulations. The MC samples are re-weighted to reproduce the distribution 
of the mean number of interactions per bunch crossing observed in the data. 


4 Event reconstruction and selection 

The following event selection criteria are applied to the events collected with a single-lepton or dilepton 
trigger. The transverse momentum and transverse energy thresholds of the single-muon and single¬ 
electron triggers are 24 GeV. Two dimuon triggers are used, one with symmetric thresholds at 13 GeV 
and the other with asymmetric thresholds at 18 GeV and 8 GeV. Eor the dielectron trigger the symmetric 
thresholds are 12 GeV. Eurthermore, there is an electron-muon trigger of thresholds at 12 GeV (electron) 
and 8 GeV (muon). 
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A primary vertex reconstructed from at least three tracks, each with > 0.4 GeV, is required. For events 
with more than one primary vertex, the vertex with the largest X p\ of the associated tracks is selected. 

Electron candidates are reconstructed from a combination of a cluster of energy deposits in the elec¬ 
tromagnetic calorimeter and a track in the ID. They are required to have pj > 1 GeV and I 77 I < 2.47. 
Candidate electrons must satisfy a loose set of identification criteria based on a likelihood built from 
parameters characterising the shower shape and track association as described in Ref. [32]. 

Muon identification is performed according to several criteria based on the information from the ID, the 
MS, and the calorimeter. The different types of reconstructed muons are: a) Combined (CB), which is 
the combination of tracks reconstructed independently in the ID and MS; b) Stand-Alone (SA), where 
the muon trajectory is reconstructed only in the MS; c) Segment-tagged (ST), where a track in the ID is 
associated with at least one local track segment in the MS; and d) Calorimeter-tagged (CaloTag), where a 
track in the ID is identified as a muon if it is associated with a minimum ionizing particle’s energy deposit 
in the calorimeter. 

The acceptance for both the CB and ST muons is \t]\ < 2.5, while the SA muons are used to extend the \t]\ 
acceptance to include the region from 2.5 to 2.7, which is not covered by the ID. CaloTag muons are used 
in the rapidity range I 77 I <0.1 where there is incomplete MS coverage. All muon candidates are required 
to have pj > 6 GeV. 

In order to reject electrons and muons from hadron decays, only isolated leptons are selected. Two 
isolation requirements are used, one for the ID and one for the calorimeter. For the ID, the requirement is 
that the scalar sum of the transverse momenta, X PT^ of all tracks inside a cone of AR = ^J(Ar])'^ + (A 0)2 - 
0.2 around the lepton, excluding the lepton itself, be less than 15% of the lepton pj. For the calorimeter, 
the Yj Et deposited inside a cone of AR = 0.2 around the lepton, excluding the lepton itself and corrected 
for contributions from pile-up and, in the case of electrons, shower leakage, is required to be less than 
30% of the muon pj (15% for SA muons) and 20% of the electron Ej. 

At the closest approach of a track to the primary vertex, the ratio of the transverse impact parameter do to 
its uncertainty, the do significance, must be smaller than 3.5 (6.5) for muons (electrons) to further reject 
leptons from heavy-flavour decays. The longitudinal impact parameter, |zol, must be less than 10 mm for 
electrons as well as muons (no vertex requirements are applied to SA muons). 

Selection of lepton quadruplets is done separately in each of the channels 4p,2e2p,4e, keeping only a 
single quadruplet per channel. Candidate quadruplets are formed by selecting two opposite-sign, same- 
flavour lepton pairs The two leading-px leptons of the quadruplet must have px > 20 and 15 GeV, 

respectively, while the third lepton must have pj > 10 ( 8 ) GeV if it is an electron (muon). The four leptons 
of a quadruplet are required to be separated from each other by A/? >0.1 (0.2) for same (different) flavour. 
At most one SA or a CaloTag muon is allowed in each quadruplet. The inclusion of final-state radiation 
to charged leptons follows the same approach as described in Ref. [33]. Each event is required to have the 
triggering lepton(s) matched to one or two of the selected leptons. All the selected 4{ events must lie in 
the 80 < < 1 000 GeV range. 

For each channel, the lepton pair with the mass closest to the Z-boson mass is selected as the leading 
dilepton pair and its invariant mass, mu, is required to be between 50 and 120 GeV. The sub-leading 
pair with the largest invariant mass, m 34 , among the remaining possible pairs, is selected in the invariant 
mass range 12 < m 34 < 120 GeV. In the 4e and 4p channels all possible pairs are required to have 
m{+(- > 5 GeV to reject events containing 7/i/^ ^ £'^€~ decays. The transverse momenta of the lepton 
pairs must be above 2 GeV. 
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5 Background estimation 


The dominant reducible backgrounds for this analysis are from Z +jets and tt processes and are estimated 
from data. Contributions from ZW, Zy, tZ as well as from the irreducible backgrounds from ttZ, VW, ZH 
and double-DY processes are estimated from simulation. 

The Z+jets and tt backgrounds are estimated using two different final states in data: + jjjj. and + ee, 
where {€ - e,iu) is the leading-lepton pair. The + jjjj. background arises from Z -i- jets and tt 
processes where the Z -i- jets contribution involves the associated production of a Z boson and heavy- 
flavour hadrons, which decay semileptonically, and a component arising from Z -i- light-flavour jets with 
subsequent njK in-flight decays. The background for + ee final states arises from associated production 
of a Z boson with other objects namely jets misidentified as electrons, which can be light-flavour hadrons 
misidentified as electrons, photon conversions reconstructed as electrons, or electrons from semileptonic 
decays of heavy-flavour hadrons. 

For both the (i+jiiJ. and €i+ee cases, the numbers of background events are estimated from a fit performed 
simultaneously to three mutually exclusive control regions, each of them providing information on one 
or more background components. The fit is based on the mass of the leading dilepton, m\ 2 , which peaks 
at the Z mass for the Z -i- jets component and has a broad distribution for the tt component. The three 
control regions are fit simultaneously to extract the different components of the reducible background, 
using a profile likelihood approach where the input template shapes for Z -i- jets and tt are obtained from 
simulation. The fitted yields in the control regions are extrapolated to the signal region using efficiencies, 
referred fo as fransfer factors, obfained from simulafion. Independenf validafion regions are used fo check 
fhe exfrapolafions. 

The fhree confrol regions for + fx/u background are defined based on fhe impacf paramefer significance 
and isolafion variables of fhe sub-leading muon pair and are consfrucfed as follows: 

• A heavy-flavour-enriched confrol region where af leasf one of fhe muons in fhe second pair fails fhe 
impacf paramefer significance requiremenf while fhe isolafion requiremenf is relaxed; 

• A lighf-llavour-enriched confrol region where af leasf one of fhe muons in fhe second pair fails fhe 
isolafion requiremenf buf passes fhe impacf paramefer significance cuf; 

• A tt-enriched region where fhe leading lepfon pair is made of opposife-sign and differenl-flavour 
leptons. For fhe muons of fhe second pair fhere is no charge requiremenf, fhe isolafion cuf is relaxed 
and fhe muons musf nol salisfy fhe impacf paramefer requiremenf. 

A validafion region fo check fhe ££ + mx background exfrapolafion is populafed by bofh Z -i- jefs and tt. 
The leading lepton pair is required fo fulfil fhe full selecfion criteria, while fhere is neifher isolafion nor 
impacf paramefer requiremenfs on fhe sub-leading muon pair. This region is used fo check fhe fif resulfs 
and verify fhaf fhe dafa and MC simulafion agree. 

The fhree confrol regions for £€ + ee background are defined based on fhe impacf paramefer significance, 
isolafion and elecfron idenfificafion requiremenfs on fhe second elecfron pair. In all confrol regions af 
leasf one of fhe elecfrons in fhe second pair musf nof safisfy fhe idenfificafion criteria. These regions are 
consfrucfed as follows: 

• A Z -I- jefs-enriched confrol region where af leasf one of fhe elecfrons of fhe second pair fails fhe 
frack isolafion and no calorimeter isolafion is required; 
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• An additional Z+jets-enriched control region where no charge requirement is made on the electrons 
of the second pair, while at least one of these electrons fails the impact parameter selection and no 
calorimeter or track isolation is required; 

• A tt-enriched region, where the leading lepton pair is selected from opposite-sign and different- 
flavour leptons. There is no charge requirement for the sub-leading electron pair. At least one 
of the electrons of the second pair fails the calorimeter isolation requirement and neither track 
isolation nor impact parameter requirements are applied. 

A validation region to check the €£ + ee background extrapolation is defined by removing the calorimeter 
isolation and requiring that at least one electron in the sub-leading pair fails the electron identification. 
Each candidate in the pair is required to pass the impact parameter and the track isolation selections, this 
region is used to check the fit outcome and verify that the data and MC simulation agree. 

The residual contributions from ZZ and ZW production in all control regions are estimated from simula¬ 
tion. The purity of the Z+jets and tf backgrounds in the control regions is above 95%. 

In the validation regions, the post-fit MC predictions agree with the data within the statistical uncer¬ 
tainty. 

The major uncertainties for the fitted reducible background come from the number of events in the control 
regions followed by the systematic uncertainty in the transfer factors. The latter is evaluated from the 
difference in the selection efficiency determined in data and simulation in dedicated control regions using 
leptons accompanying Z —» candidates, where the leptons composing the Z-boson candidate are 

required to satisfy isolation and impact parameter criteria. Events with four leptons are excluded. Eor 
the MC estimated background the systematic uncertainties mainly come from theoretical cross-section 
uncertainties for different processes and from luminosity uncertainties in normalisations. The differential 
distributions for all background processes are taken from simulation. 

Table 1: Numbers of expected background events for different processes and channels. 


Process 

4e 

4/1 

leljd 

tt 

0.45 ± 0.24 

0.68 + 0.19 

1.3+ 0.5 

Z + jets 

0.6 ± 0.29 

5.3 + 1.5 

6.3+ 1.4 

Diboson 

1.25 + 0.18 

0.83 + 0.18 

2.84 + 0.34 

Triboson 

0.67 + 0.12 

0.97 + 0.14 

1.46 + 0.19 

Z-i-top 

0.62 + 0.15 

1.19 + 0.32 

1.7+ 0.5 


The total number of background events estimated from data and MC simulation is 26.2±3.6. Numbers 
of background events expected per channel estimated for different processes are shown in Table 1 . The 
background estimation was cross-checked with an alternative method, described in Refs. [1, 34], called 
the fake-factor method. The results from this cross-check are found to be consistent within uncertainties 
with those described above. 


6 Cross-section extraction method 

Two cross sections are extracted from the number of observed events. One is the fiducial cross section, 
cr^^, in the experimental phase space defined by the event selection criteria and the other is the cross 





section, cr“^, in an extended common phase space where electrons and muons have the same geometric 
and kinematic acceptance. The fiducial phase space is defined in Table 2. The exfended phase space for 
fhe cross-seclion exfracfion is defined by 80 < niin < 1 000 GeV, m^+i- > 4 GeV, > 2 GeV, and 
fhe presence of four leptons each wifh pj > 5 GeV and \r]\ < 2.8. 

The cross section measuremenf is performed using a likelihood til described below. For a given channel i, 
the observed number of events, follows a Poisson distribution, PoisCA^Q;,^, the mean of which, 

Vpred “ + A^b> is ih® sum of the expectations for signal and background yields. These yields depend on 

the fiducial cross section and fhe nuisance paramefers, x, which represenf fhe experimenfal and fheorefical 
uncerfainfies as: 


= Ni(crlf,0){\+Y,XkS[), ( 1 ) 

k 

Nl(£) = V'(0)(l + 2]x,4), (2) 

k 

where S ^ and are fhe relative sysfemafic effecfs on fhe signal and background, respectively, due fo fhe 
k-fh source of sysfemafic uncerfainfy. The cenfral expecfafion of fhe signal yield, corresponding fo fhe 
sysfemafic sources af fhe nominal value (referred fo as fhe nuisance-free expecfafion), is given by : 

K{(Tl^,0) = £-C^k-Kr-all (3) 

where £ is fhe infegrafed luminosify, and is fhe rafio of fhe number of accepfed signal evenfs fo fhe 
number of generated evenfs in fhe fiducial phase space. Corrections are applied fo C\e fo accounf for 
measured differences in frigger and reconsfrucfion efficiencies befween simulafed and dafa samples. The 
Cu values are 53.3%, 82.2% and 67.7% for fhe 4e, 4p, and 2e2p channels, respectively. The confribufion 
from T-lepfon decays is accounfed for by a correcfion term K-^ - I , where is fhe number 

of accepfed simulafed 4£ evenfs in which af leasf one of fhe Z bosons decays info r-lepfon pairs, and 

Slg 

is fhe number of accepfed simulafed ZZ evenfs wifh decays info elecfrons or muons. 

Cross-section measuremenfs are exfracfed for a single channel or any combinafion of channels, using a 
likelihood mefhod. The likelihood function is: 

L(<, f) = f[ Pois(V',„A)) • e-^, (4) 

i 

where fhe producf runs over fhe channels fo be considered. 

For fhe exfended phase space fhe likelihood function is parameferised as a function of fhe exfended cross 
section similar to fhe one shown in Eq. (3) and multiplied by fhe fiducial accepfance A^e, which is fhe ratio 
of fhe number of evenfs wifhin fhe fiducial phase-space region fo fhe fofal number of generafed evenfs in 
fhe exfended phase space. The fiducial accepfance is evaluated using simulation fo be 41.6%, 50.3%, 
and 42.2%, for fhe 4e, 4p, and 2e2p channels, respectively. The differences are due fo fhe elecfron and 
muon geomefric defecfion coverage. 

To find fhe cenfral value of fhe cross secfion cr, the likelihood function is maximised simultaneously with 
respect to the nuisance parameters and cr. Correlations between the signal and background systematic 
uncertainties are taken into account in the likelihood fitting procedure. 
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Table 2; List of selection cuts which define the fiducial region of the cross-section measurement. Same-flavour 
opposite-sign lepton pairs are denoted as SFOS, the leading lepton pair mass as m\ 2 , and the sub-leading lepton 
pair mass as ^ 34 . The four-momenta of all final-state photons within AR = 0.1 of a lepton are added to the four- 
momentum of that lepton. 


Lepton selection 

Muons: 

Pt > 6 GeV, \r]\ < 2.7 

Electrons: 

Pt > V GeV, IpI < 2.5 

Lepton pairing 

Leading pair: 

SFOS lepton pair with 
smallest \mz - ma\ 

Sub-leading pair: 

The remaining SFOS 
with the largest 

For both pairs: 

pf > 2 GeV 

Event selection 

Lepton 

> 20,15,10(8 if p) GeV 

Mass requirements: 

50 < mi 2 < 120 GeV 

12 < m 34 < 120 GeV 

Lepton separation: 

AR{€i,€j) > 0.1 (0.2) 
for same- (different-) 
flavour leptons 

J/ifj veto: 

> 5 GeV 

A£ mass range: 

80 < m 4 i < 1 000 GeV 


7 Systematic uncertainties 


Systematic uncertainties on the measurement arise from uncertainties on the integrated luminosity, the 
experimental calibrations of the lepton energy and momentum, and the lepton detection efficiencies, as 
well as the theoretical modelling of signal acceptance, and the background estimation. The overall un¬ 
certainty on the integrated luminosity is ±2.8%, which is derived following the same methodology as 
that detailed in Ref. [35]. A summary of the relative uncertainties of and A 4 ^ x is given in 

Tables 3 and 4. 

The effect on the expected signal event yields due to experimental systematic uncertainties is determined 
from the uncertainties on lepton energy and momentum scales and resolutions, as well as the uncertainties 
on efficiencies of the lepton reconstruction and identification. The major contributions come from the 
uncertainties on lepton reconstruction and identification efficiencies [36-38]. 

The uncertainties on the signal acceptance for both Cu and A 4 ^ include theoretical uncertainties from the 
choice of QCD scales and PDF set. The scales are varied independently from 0.5 to 2.0 times the nominal 
values of /ir and //p. The PDFs uncertainties are estimated by using the envelope [39] of variations of 
different PDF sets, CTIO, MSTW2008 [40] and NNPDF2.3 [41]. 

The Cd,t uncertainty is mostly experimental and of the order of 2-5%, while the A 4 ; uncertainty is entirely 
theoretical and of the order of 3-5%. A range of values of the relative uncertainties on the are 
given by 4.9%, 1.9%, and 2.5% for the Ae, Ajj., and 2e2/i, respectively. The uncertainties on €41 due to 
higher-order corrections to the gg production processes are less than 0.6%. This is estimated by applying 
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an approximate NNLO ^-factor determined for the Higgs-boson production [23], assuming that it is 
applicable to the normalisation of the continuum gg ZZ production cross section. 

Uncertainties on C^t, as a function of mm and are also computed for the differential cross section 
measurements. In the mass region {mm < 150 GeV), the relative uncertainties on Cm vary in the range 
of 4-9%, 1.7-2.7%, and 2-5% for the \e, 4p, and 2e2p channel, respectively. In the mass region mm > 
150 GeV, they are almost constant as a function of mm and are about 4%, 1.8%, and 3% for the 4e, 4p, 
and 2e2p channel, respectively. 

The relative uncertainties on Am are 1.2%, 1.0%, and 1.6% for the 4e, 4p, and 2e2p channel, respectively, 
evaluated by comparing POWHEG-BOX and MCFM MC samples with the same approach for the QCD 
scales and the PDF uncertainties as described earlier. The QCD scale uncertainties do not change when 
going from NFO to NNFO for the signal normalisation for the qq —> events [16]. An additional 
uncertainty (3-4%) is included in the Am uncertainty estimate to account for the uncertainty of the Higgs- 
boson NNFO X-factor normalisation correction of the non-resonant signal from gluon fusion (labelled 
“extra gg corrections” in Tables 3 and 4). 


Table 3: The combined relative uncertainties on the efficiency correction factor Cm, evaluated as the sum in quad¬ 
rature of the uncertainties from different sources, including electton and muon identification and theoretical uncer¬ 
tainties due to PDFs, QCD scales, and parton shower modelling. Extra uncertainties due to higher-order corrections 
for the gg process (NNLO X-factors for Higgs-boson production applied to the inclusive gg process) are also given. 


Sources 

ACmICm 

Ae 4jU 2e2// 

Experimental (e) 

4.8% 

- 

2.3% 

Experimental (jj.) 

- 

1 .8% 

0.9% 

Theoretical 

0 .1% 

0 .1% 

0 .2% 

Extra gg corrections 

0 .6% 

0 .2% 

0.3% 

Combined uncertainty 

4.9% 

1.9% 

2.5% 


Table 4: Theoretical relative uncertainties on the fiducial acceptance Am and Am x Cm due to PDFs, QCD scales and 
parton shower modelling. Extra uncertainties due to higher-order corrections for the gg process (NNLO X-factors 
for Higgs-boson production applied to the inclusive gg process) are also given. 


Sources 

4e 

AAml^u 

An 

le2n 

Theoretical 

1 .2% 

1 .0% 

1 .6% 

Extra gg corrections 

4.0% 

3.0% 

3.9% 

A{Am X Cm)!{Am x Cm) 

Theoretical 

1.4% 

1 .1% 

1.7% 

Extra gg corrections 

4.6% 

3.2% 

4.2% 


The overall uncertainty on the background estimation is ±14%. The contributions from different sources 
and channels are given in Table 1. 
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Table 5: Summary of the observed and predicted numbers of events in different Ai channels. denotes the 
selected number of data candidates. denotes the total predicted number of events (including t contributions) 

for signal plus background. and the denote the predicted non-gg signal and the gg signal (no NNLO 

A'-factor has been applied), respectively. A^^*- denotes the t contributions. A^bkg denotes the total estimated number 
of background events (from data and MC simulation). The listed uncertainties of the expected number of signal 
events include statistical and experimental systematic uncertainties. 


Channel 

^Data 

^Total 

'^expected 

^signal 

^signal 

^MC 

f^bkg 

Ae 

85 

80 + 4 

68.4 + 3.4 

6.24 + 0.31 

1.28 + 0.06 

3.6+ 0.5 

Ap 

156 

150.2 ± 2.9 

128.2 + 2.5 

11.00 + 0.21 

2.18 + 0.09 

9.0+ 1.5 

lelp 

235 

205 + 5 

172 + 5 

16.0 + 0.4 

3.08 + 0.13 

13.6 + 2.1 

Tofal 

476 

435 + 9 

369 + 9 

33.3 + 0.8 

6.54 + 0.14 

26.2 + 3.6 


8 Results 

8.1 Cross-section measurements 

The numbers of expected and observed events after applying all selection criteria are shown in Table 5. 
A total of 476 candidate events is observed with a background expectation of 26.2 ± 3.6 events. The 
observed and predicted and p‘^ distributions for the selected events are shown in Figure 3. 

The measured cross sections in the fiducial and extended phase space for different At channels are sum¬ 
marised in Table 6 and compared to the SM predicted cross sections. The combined At cross section in 
the extended phase space is found to be 73 +4 (stat.) ±4 (syst.) +2 (lumi.) fb, compared to a SM prediction 
of 65 ± 4 fb. One should note that the cross section for non-resonant ZZ production from the g'^-induced 
signal is only calculated at LO approximation, which could be significantly underestimated. 


Table 6: Measured cross sections in the fiducial phase space (cr®‘*) and extended phase space (cr“‘), compared to 
their SM predictions (calculations described in Section 3). One should note that the non-resonant gg-induced signal 
cross section is only calculated at LO approximation. 


At 

Measured cr^‘^ [fb] 

SM cr“ [fb] 

Measured cr®’“ [fb] 

SM 0 -“' [fb] 

\e 

7 . 4 :°;® (stat) (syst) (lumi) 

6.9 ± 0.4 

17.8];|(') (stat) ];[ j (syst) (lumi) 

16.4+1.0 

Ap 

8.7];q 7 (stat) (syst) (lumi) 

8.3 ±0.5 

17.31[;^ (stat) (syst) (lumi) 

16.4 ± 1.0 

lelp 

15.9l[;[ (stat) (syst) (lumi) 

13.7 + 0.9 

37.7+1^ (stat) +1^ (syst) +[;[ (lumi) 

32.1 + 2.0 

Total 



73];4 (stat) (syst) ^2 (lumi) 

65 + 4 


8.2 Differential cross-section measurement 

The measurement of the differential cross-section is performed in the fiducial phase space defined in 
Table 2. The evenfs from all fhree At channels are combined info a common sample for fhe unfolding 
procedure. The unfolding is done as a funcfion of fhe fwo kinematic variables and . The 014 ^ 
specfrum is essential for fhe sfudy of fhe differenf producfion mechanisms, while fhe pt^ specfrum is 
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sensitive to higher-order QCD corrections and to QCD resummation effects at small p‘^ [10]. The high- 
region is sensitive to top-loop effects in gg ^ H production as well as to anomalous triple-boson 
couplings. 

The iterative Bayesian unfolding [42] is applied here. In the unfolding of binned data, the effects of the 
experimental acceptance and resolution are expressed in terms of a response matrix, where each element 
corresponds to the probability of an event in the /-th generator level bin being reconstructed in the y'-th 
measurement bin. The response matrix is combined with the measured spectrum to form a likelihood, 
which is then multiplied by a prior distribution to produce the posterior probability of the true spectrum. 
The SM prediction is used as the initial prior, and once the posterior probability is obtained, it is used as 
the prior for the next iteration. The spectrum becomes insensitive to the initial prior after a few iterations. 
The differences befween successive iferafions are used fo esfimafe fhe sfabilify of fhe unfolding mefhod. 
In fhis analysis four iferafions are performed. 

The unfolded disfribufions are shown in Figure 4, where fhe differential cross secfion is presenfed as a 
funclion of m 4 f and p^^ and compared fo fheory predicfions. The dafa points shown in the figures are 
fhe measurements with combined statistical and systematic uncertainties. The theoretical predictions 
are the sum of the differential cross sections of the qq and gg A€ processes. The LO cross 

sections are used for the non-resonant ^^-induced signals. The cross sections of the on-shell Higgs boson 
are normalised to include the NNLO QCD and NLO EW effects as summarised in Ref. [19]. The qq 
production processes are corrected with NNLO QCD and the NLO EW /T-factors for the m^t spectrum 
for > 2 X mz- Lor the pt^ spectrum, the qq signal prediction is calculated by POWHEG-BOX at 
NLO. 

The uncertainties on the differential cross-section measurements are dominated by the statistical uncer¬ 
tainties of the data. Lor example, in the regions between the Z and Higgs boson peaks and between 
the Higgs-boson mass mu and = 180 GeV, the statistical uncertainties are of the order of 45% and 
20%, respectively. In the high-mass region {m^i > 180 GeV) they are of the order of 10%. Lurthermore, 
one should note that the NNLO QCD corrections are not available for the qq production calculation 
for the mass region < 2 x niz- 

In the bin of 120-130 GeV, which is dominated by the resonant Higgs-boson contribution, the ratio 
of data to the MC prediction is compatible with the ATLAS measurement [33] of the Higgs-boson signal 
strength of pn - The data points in the spectrum between 140 to 180 GeV are slightly 

more than Icr above the theoretical predictions, where the NNLO QCD correction is not yet available. 
Some discrepancy is also observed in the lowest bin and in the region between 30-50 GeV of the p‘}^ 
spectrum. 


8.3 Extraction of the gg signal contribution in the m^t > 180 GeV region 

The extraction of the signal strength of the non-resonant gg Ai production is performed in the high- 
mass region, >180 GeV, where this production mode is dominated by the continuum gg ZZ 
process through a quark-box diagram intermediate state (see Ligure 1(d)). Additional contributions come 
from the off-shell Higgs-boson production and the interference between Higgs boson and continuum ZZ 
production. 

The m^f: spectrum is chosen as the discriminant to extract the gg signal strength with respect to the LO gg 
prediction: pgg = o-(data)/cr(LO). 
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The contribution of the qq ZZ production is constrained to the best theory knowledge (which accounts 
for QCD NNLO and EW NLO m 4 /-dependent corrections) and jigg is extracted from a likelihood fit using 
the reconstructed distributions. The experimental uncertainties are treated as fully correlated between 
qq and gg processes. The theoretical uncertainties, including the uncertainties on the normalisation of the 
qq ZZ —> Tf”, the shapes of spectra from both the qq and gg initial states, and the acceptance, are 
taken into account. The distribution of the data, the fit, the expectation from non-^^ signal processes 
and the background are shown in Figure 5. The fit result is= 2.4± 1.0(stat.)±0.5(syst.)±0.8(theory). 
This result corresponds to a ^g-initiated cross section of 3.1 fb, which has the same relative uncertainties 
as Hgg itself in the inclusive fiducial volume as defined in Table 2 wifh fhe addifional requiremenf of 
m 4 ( >180 GeV. The largesl uncerfainfy is sfafisfical. The fheorefical uncertainly is mainly due fo fhe 
normalisalion uncertainly of fhe qq —> ZZ process. 

The fheorefical eslimale of m 4 ^-dependenf /f-faclor for off-shell Higgs boson produclion given in Ref. [23] 
is in a range of 2.7-3.1 (wifh CTIONNLO PDF) and lhal given in Ref. [24] for fhe interference lerm is 
2.05-2.45. These fheorefical sludies confirm lhaf fhe gluon sofl-collinear approximafion predicls similar 
/f-faclors for off-shell Higgs-boson and inferference, hence supporfing fhe assumpfion of a similar K- 
facfor for fhe confinuum ZZ produclion. These fheorefical calculafed A'-facfors are compalible wifh fhe 
resull obfained by Ihis analysis, where fhe g'^-iniliated A€ evenls are produced predominanfly from fhe 
confinuum ZZ produclion. 

Applying fhe higher-order corrections fo bofh fhe cross section of fhe off-shell Higgs-boson production 
and fhe conlribulion of fhe inferference lerm, while keeping fhe FO cross section for fhe confinuum 
gg ZZ produclion, fhe change of fhe Hgg til resull is negligible (approximately, = 0.01). 


14 



Figure 3: Data and MC prediction comparison for selected events as a function of the invariant mass m 4 { (top) and 
the transverse momentum (bottom) of the four-lepton system. The solid colours show the expected contributions 
from signal and background and the black points represent data with statistical error bars. 
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Figure 4; The measured differential cross-section distributions (the black points) of(left) and (right), unfol¬ 
ded into the fiducial phase space, and compared to theory predictions (red histogram). The combined statistical and 
systematic uncertainties of the measurements are shown as the error bars of the unfolded spectra. The theoretical 
predictions are the sum of the differential cross sections of the qq —> and gg —> processes, where the LO 
cross sections are used for the non-resonant grgr-induced signals, and the cross sections of the on-shell Higgs boson 
and the qq production processes are corrected with the NNLO /T-factors for the m/ii spectrum; except for the 
where only the NLO and LO predictions are used for the qq and the gg processes, respectively. The total theoretical 
uncertainties are shown as error bands evaluated by the sum in quadrature of the contributions from parton showers, 
QCD scales, PDF sets, and electroweak corrections. 
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Figure 5: Comparison of the m 4 t spectra between the data (black points with error bars) and the prediction (red 
histogram) after the likelihood fit of jigg. The non-gg signal from the theoretical prediction (blue histogram) and the 
background (brown histogram) are also shown. The gg contribution is the difference between data and the sum of 
the non-gg signal and the background. 
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9 Conclusion 


The measurement of four-lepton produetion in proton-proton eollisions at = 8 TeV is presented using 
data eorresponding to an integrated luminosity of 20.3 fb“^ eolleeted with the ATLAS deteetor at the 
LHC. In total, 476 eandidate events are observed, with a baekground expeetation of 26.2 ± 3.6 events, 
in the four-lepton invariant mass range between 80 and 1 000 GeV. The M produetion eross seetions 
are determined in both fidueial and extended phase spaees. The measured eross seetion in the extended 
phase spaee, defined by 80 < mm < 1 000 GeV, m^+i- > 4 GeV, > 2 GeV, four leptons eaeh with 
Pt > 5 GeV and \r]\ < 2.8, is found to be 73 ± 4 (stat.) ±4 (syst.) ±2 (lumi.) fb, and is eompared to a SM 
predietion of 65 ± 4 fb. The measurements of the 4£ differential eross seetions are performed by unfolding 
the mm and the speetra. In the mass range above 180 GeV, assuming the theoretieal eonstraint on the 
qq produetion eross seetion ealeulated with perturbative NNLO QCD and NLO eleetroweak eorreetions, 
the signal strength of the gluon-fusion eomponent with respeet to the LO predietion is determined to be 
Pgg = 2.4 ±1.0 (stat.) ± 0.5 (syst.) ± 0.8 (theory). 
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